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Abstract 

OS 

We study the multi-wavelength signal induced by pairs annihilations at the Galactic 
center (GC) of a recently proposed dark matter (DM) candidate. The weakly in- 
teracting massive particle (WIMP) candidate, named ^4_, is the first Kaluza-Klein 
O mode of a five dimensional Abelian gauge boson. Electroweak precision tests and 

the DM cosmological bound constrain its mass and pair annihilation rate in small 
ranges, leading to precise predictions of indirect signals from what concerns the par- 
ticle physics side. The related multi-wavelength emission is expected to be faint, 
unless a significant enhancement of the DM density is present at the GC. We find 
that in this case, and depending on few additional assumptions, the next generation 
of gamma-ray and wide-field radio observations can test the model, possibly even 
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O ' with the detection of the induced monochromatic gamma-ray emission. 
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In the last decades gravitational evidences for dark matter (DM), based on dif- 
ferent observables, such as rotation curves, velocity dispersions, gravitational 
lensing, large scale structure maps and cosmic microwave background (CMB) 
anisotropics have been accumulated at the galactic, cluster and cosmological 
scales. Still, the identification of the DM component remains one of the most 
challenging issue of the physics today (for recent reviews see |H;0)). Weakly 
interacting massive particles (WIMPs) are a well motivated class of candi- 
dates for the non-baryonic component. The WIMP paradigm is well-known: 
In thermal equilibrium in the primordial bath, WIMPs decouple in the non- 
relativistic regime and the weak interaction leads the relic abundance to be 
of the order of the mean energy density of DM in the Universe today. Being 
(weak) interacting particles, WIMPs can annihilate in pairs in astrophysical 
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Fig. 1. Left Panel: For a few selected values of the DM mass, isolevel curves of the 
parameter space (A& = (M&_ — Ma_)/Ma_, p s ) in which the A_ relic abundance 
matches the best fit value fiflM^ 2 from cosmological observations. The solid and 
dashed lines refer to the two different setups described in the text. Right Panel: 
Gamma-ray and electron/positron differential spectra per annihilation of the DM 
candidate A- in the minimal framework. 



structures, inducing detectable signatures, such as antimatter and neutrino 
fluxes, and multi-wavelength spectra. Complementary to direct DM searches 
and to collider experiments testing extensions to the standard model (SM) of 
particle physics embedding a WIMP candidate, indirect detection can provide 
crucial information about the fundamental nature of DM. 



In Ref. (|3|) a new viable DM candidate was introduced. It is embedded in a flat 
five dimensional (5D) model compactified on an orbifold 5 1 /Z 2 , withgauge- 
Higgs unification and explicit Lorentz symmetry breaking in the bulk (|4|; la; |6|) . 
The DM candidate, named A_, is the first Kaluza-Klein (KK) mode of an 
Abelian 5D antiperiodic gauge boson field. It is the lightest KK particle odd 
under the "mirror symmetry", a discrete Z2 symmetry introduced in Ref. (0) to 
improve the naturalness of the model, namely to reduce the fine-tuning needed 
to stabilize the electroweak (EW) scale. In a relevant fraction of the parameter 
space, the DM candidate is nearly degenerate with the lightest antiperiodic 
fermion 6_ and the antiperiodic gluon g_ (see the mass spectrum in Fig. 1 
of (0)). The mass for the first KK mode of antiperiodic fields is given at tree 
level by: 
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where /_ denotes an antiperiodic fermion, R is the compactification radius of 
the covering circle S , pi and kf are the Lorentz breaking parameters of gauge 
boson and fermion respectively, and my is the fermionic bulk mass. EW bounds 
force the compactification scale of the model in the Multi-TeV regime (0) , lead- 
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ing to a lower bound for the DM mass: Ma_ > 2.35 TeV, assuming p ~ 1. The 
A- pair annihilation rate is quite small (cr ann t> < 5- 10~ 28 cm 3 s _1 ) compared to 
WIMPs in more standard scenarios (e.g. the lightest neutralino in supersym- 
metry and the B 1 in universal extra dimension (UED) scenario, see Ref. (Ill)); 
however, being 6_ and g_ strongly interacting particles, coannihilation effects 
greatly enhance the effective annihilation rate, leading to relic abundances 
allowed by cosmological observations. The curves in Fig. [Tk show the mass 
splitting between the DM candidate and the coannihilating particles, such 
that the A_ relic density matches the cosmological amount of DM today (0). 
They are expressed in terms of (A b , p s ) with A b = {M b _ — Ma_)/Ma_, tak- 
ing into account radiative corrections as described in Ref. (0). Fig. [Ih shows 
two different setups. Indeed, the gauge group of the model is G x G', where 
G = SU(3) W x SU(3) S x U(l) t and G' = SU(3)' S x U(l)' as in the first con- 
struction of Ref. pa) (dashed lines) or G' = U(l)' as in the framework (in 
the following called "minimal" framework) considered in Ref. (0) (solid lines); 
antiperiodic gluons are present only in the first setup and differences between 
the two cases for what concerns the dark matter relic density computation are 
described in Ref. Combining EW and cosmological bounds, the A_ mass 
is constrained in the narrow window 2.35 - 5 TeV. Naturalness arguments can 
restrict even more the parameter space of the model. Indeed the value for the 
mass preferred by EW observables is ~ 3 TeV (@) and the fine tuning asso- 
ciated to the DM relic density is minimized by the minimal framework (0), 
where A_ annihilates only with 6_, leading to (M b _ — Ma_)/M a _ < 7% (see 
Fig. Eh). 

Couplings with SM fermions are highly suppressed^ since the latter (with the 
exception of bottom and top quarks) are mainly localized on the 4D brane at 
y = 0, where the A_ wavefunction vanishes, being antiperiodic on S . As we 
can see from the mass spectrum in Ref. j3) , some non-standard states are en- 
ergetically accessible by the A_ pair annihilation. They belong to SU(2) W sin- 
glet, doublet or triplet representation and they are the first massive eigenstates 
obtained from the diagonalization of the Lagrangian describing the coupling 
between 5D periodic fermions and 4D localized fermions. Non-SM singlets 
and doublets are the KK excitations of the SM fermions and they are mainly 
localized. The case of triplets is different: they are not coupled with boundary 
fields and their wavefunctions are constant in the bulk, hence largely over- 
lapping the A_ wavefunction. Moreover the associated number of degrees of 
freedom is huge and so the dominant final states of the A_ annihilation cross 
section are triplets fermions. More precisely the annihilation branching ratios 
are: 75% into b + b + , 24% into r + f + (with b + and r + being the SU(3) W triplet 
of the bottom and tau tower, respectively) and 1% into SM quarks. The sub- 



1 This fact implies a very small elastic scattering cross section between A_ and light 
quarks, and the expected direct DM signals are well below the sensitivity of current 
detection experiments. 
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sequent decays of b + and r + generate quark pairs (38%), r lepton pairs (6%) 
and neutrinos (6%), charged (25%) and neutral (12%) weak gauge bosons, and 
Higgs bosons (12%). 

Photons cannot be directly produced by WIMP pair annihilations at tree level. 
A continuum photon spectrum is generated in cascades with hadronization of 
unstable two-body annihilation final states into n°s and their subsequent de- 
cays. The energy of these photons is in the 7-ray band. Electrons and positrons 
can be directly or indirectly produced in WIMP pair annihilations as well. The 
associated radiative processes can act as source for a multi-wavelength spec- 
trum covering from radio to soft-gamma ray frequencies. In Fig. [TJd we show 
the differential energy spectra per ^.-annihilation into 7-rays and e + — e~ in 
the minimal DM framework (variations of the DM mass within the allowed 
range do not affect the spectra in a sizable way). In the first case, on top of the 
spectrum originated from ir° decay, we consider the contribution of primary 
gamma-rays from final state radiation following the line of Refs. (Q) and (0). 
We derive the differential yields through simulations of decay and hadroniza- 
tion performed with the PYTHIA Monte-Carlo code (fiol ). The two spectra 
are soft since quarks and gauge bosons are the dominant annihilation modes. 
From the point of view of indirect searches, this feature distinguishes A_ from 



the UED WIMP candidate B 1 dill; ll2|) , whose pairs annihilation branching 



ratios are dominated by charged leptons and harder spectra are produced. The 
electron/positron and gamma-ray yields of Fig. [TJd are at a comparable level, 
being mostly generated by the production and decays of charged and neutral 
pions, respectively, with the two chains having analogous efficiency. 

As listed in Table 2 of Ref. ([^), all the ^.-annihilation processes occur through 
t or u-channels mediated by an antiperiodic fermion. At a given DM mass, 
the only free parameter affecting in a sizable way the cross section computa- 
tion is the mass of the mediator. As already mentioned before, in the minimal 
DM framework, the relic abundance is driven mostly by the 6_ coannihila- 
tion, highly constraining M b _ and hence k b (see Eq. [TJ). It leads the total A_ 
annihilation cross section within a small range, since the triplet pairs associ- 
ated to the 6-multiplet are the dominant annihilation modes. The 5D Lorentz 
symmetry breaking was introduced to achieve the correct value for the top 
mass; the Lorentz breaking parameter ki associated to other fermions can be 
safely taken ~ 1. For our purposes k b and k T are relevant in the computation 
of the annihilation cross section in the non-minimal scenario, where we as- 
sume h < 2. The allowed total annihilation cross sections as a function of the 
WIMP mass are shown in Fig. [TJby the filled band; in the minimal framework 
this region shrinks to its upper boundary. 

The WIMP candidate we consider has a small annihilation rate and a quite 
heavy mass. Looking for a WIMP induced signal might seem hopeless, unless 
we concentrate on a region where the DM density is very large and where 
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emissions from other astrophysical sources are faint. In the following we focus 
on photon emissions at the Galactic Center (GC), where the first condition is 
definitely fulfilled, while the second is not completely satisfied, at least in the 
wide range of frequency over which a source labeled Sgr A* has been detected. 
On other hand, the luminosity associated to this source turns out to be quite 
low, at a level comparable to the expected DM signal. 



In the determination of the v4_-induced emission, the ingredients related to 
the particle physics side are quite strictly constrained, while the astrophysi- 
cal uncertainties remain large, as for any WIMP model. The Milky Way DM 
halo profile is not completely established, in particular for what concerns its 
inner region. As well known, there is some tension between N-body simula- 
tion results and observations of spiral galaxies about the central cusp/core. 
Sufficiently bright DM signals, at the level needed to test the WIMP model 
under investigation, appear only for halo profiles with a large overdensity in 
the GC region. We follow the analysis in Ref. (fl3l ) and focus our attention on 
a distributions obtained from the evolution of a Navarro-Frenk- White (NFW) 
profile (Il4l). including the deepening in the Galactic potential well generated 
by the slow adiabatic formation of the supermassive black hole (SMBH) and 
of the stellar component in the inner Galaxy. In such profile, hereafter labeled 
A sp , the effect of self annihilation triggers the density in the innermost reg ion. 
The numerical profiles were kindly provided by the authors of Ref. (1131 ) for 
values of (av) and M DM in the filled band in Fig. [U while we derive scaling 
laws in the rest of the parameter space starting from results reported again in 
Ref. (flT 



In the GC case, synchrotron radiation and, to a smaller extent, inverse Comp- 
ton scattering on CMB and starlight are the most significant radiation mecha- 
nisms. The estimation of the associated emissions requires a description of the 
electrons/positrons propagation and of their energy transfers after production 
in WIMP annihilations. We need a model of the galactic medium, to fix the 
diffusion coefficient, the magnetic field, the advective/convective velocity and 
the absorption effects. For this treatment we refer to Ref. (|15l ). 



The compact radio source Sgr A* associated to the SMBH has been detected 
at the GC (see the catalog in Ref. ( fl6l ) and reference therein), together with 
its infrared (1 1.71) and X-ray (|18l ) counterparts. A gamma-ray emission from 
the GC region has been detected as well (1 19 ; l20l ) , but with experimental an- 
gular resolutions not sufficient to identify the source and its precise location. 
Any of the portions of such multi-wavelength spectrum turns out to be in- 
compatible with emissions induced by WIMP annihilations and the detected 
signals will be exploited here to derive upper limits. Also diffuse emissions 
from the inner region of the Galaxy has been detected at different frequency 
bands, and, in case of shallow DM halo profile, can severely constrain WIMP 
models. The procedures implemented to extract the limits shown in Figs, dh 
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Fig. 2. Exclusion limits on the A- annihilation cross section as a function of the 
WIMP mass. The Left and Right Panels show the cases of A sp and NFW profiles, 
respectively. 



and db were outlined in Ref . (Il5l ) . We assume A_ accounting for the whole 
DM content of the Universe and all the numerical calculations are performed 
with the help of the DarkSUSY package (12 ll ). Together with bounds associated 



to the mostly investigated profile in this Letter, i.e. the A sp profile, we com- 
pute, for comparison, limits on the WIMP parameter space in case of a NFW 
profile, namely the mostly investigated case in the literature. We plot the 
tightest bounds in gamma-ray and radio bands obtained from spectral and 
angular analysis, comparing the WIMP signals with the emission detected 
by the 7-ray air Cherenkov Telescope (ACT) HESS (I20l ; 1221 ) and with upper 
bounds in the radio surveys of Refs. (|23l ) and (1241 ) . In the X-ray band, syn- 
chrotron emission would require very strong magnetic field, especially in case 
of soft electron/positron spectrum. This could be possible only in the inner- 
most region of the Galaxy, depending on the model considered for accretion 
flow around SMBH, hence the size of the DM induced source is very small. 
Limits on WIMP parameter space can be extracted by the comparison with 
the Sgr A* emission detected by the Chandra observatory (fl8l ). but they are 
highly model dependent. We plot the weakest constraint among the three cases 
with different choice of magnetic field radial profile of Ref. (Il5l ) . The angular 
size of the emission induced by the inverse Compton scattering on CMB is 
much larger and the signature estimate involves more reliable assumptions on 
the magnetic field strength at larger scales. The limit extracted by the com- 
parison with the detected X-ray diffuse emission ( 1251 ) (dashed-dotted lines) 
is much less constraining (but more robust) with respect to the limit asso- 
ciated to the point-like synchrotron source (dotted lines); the fact that the 
latter is excluding the whole A_ parameter space in the A sp case should not 
be overemphasized, given the critical extrapolations involved in this result. 

Then we derive projected limits from forthcoming gamma-ray surveys and 
wide-field radio observations. For heavy WIMP models, the parameter space 
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can be more efficiently studied by ACTs rather than space satellites, due to the 
different energy ranges of detection. We consider the next generation of ACT, 
under development by the Cherenkov Telescope Array (CTA) project, and 
scheduled for 2013, assuming performances as outlined in (0). A diffuse radio 
emission was reported both in the Milky Way atlas of Ref. (1271 ) and in the GC 
image of Ref. (|28l ). However, the two surveys have too poor angular resolutions 
to resolve the spatial profile of the emission in the innermost region. In the 
GC map of Ref. (0), such emission does not seem completely isotropic and 
tight constraints are derived from patches of the map with no astrophysical 
background. The radio projected limits plotted in Figs.dh anddb are extracted 
again following Ref. (Il5l ). but assuming a detector sensitivity improved by a 
factor 10, as proposed in the EVLA project (1291 ) . 



The scale at which the formation of the SMBH could have influenced the DM 
distribution is far below the resolution of both numerical simulations and ob- 
servations. The related DM spike in the A sp profile greatly enhances signals 
in the innermost region of the GC and the comparison with the Sgr A* source 
is very constraining, especially for (<tv)/Mdm > 10~ 32 cmT^ s~ x GeV~ x (Il3l ). 
The limits associated to diffuse emissions are less constraining, since involve 
angular scales where the enhancement in the DM distribution A sp is less pro- 
nounced with respect to an NFW profile, being related to the deepening in 
the potential well induced only by the stellar component. For the same reason, 
being the DM induced radio source more extended than the DM source itself, 
and thus than the gamma-ray source, the bound associated to wide field radio 
signal is less stringent with respect to gamma-ray limit in case of A sp profile. 
The picture is reversed for the NFW distribution. In case of A sp profile, all 
the multi-wavelength constraints extracted from past surveys, excluding the 
synchrotron X-ray bound, do not limit the region allowed by cosmological and 
EW bounds (filled band). On the other hand, in the next decade, the model 
could be completely tested through its gamma-ray emission by the CTA ex- 
periment. The plotted exclusion curve is computed assuming an effective area 
A e ff = 1 km 2 and an exposure time t exp = 250 hours in 5 years of collect- 
ing data. Depending on the properties of the galactic radio diffuse emission 
at small scales, the EVLA project could test the v4_ radio profile in a large 
fraction of the parameter space, covering basically the whole region of the 
minimal DM framework. In the case of NFW profile, no significant constraint 
can be derived. Note however that radio wide field observations can be much 
more efficient than gamma-ray measurements. 

Radio observations with a wide field of view have detected extended emis- 
sions from the GC region. In Fig. [3h we plot schematic representations of the 
angular shape of the signals at 90 cm, as detected in the map of Ref. (j2~j ) 
(FWHM=43") and Ref. (0) (FWHM=40'). For both we sketch the profile 
of the extended source along its longitudinal axis. The level of the DM in- 
duced emission filtered over the same experimental angular resolutions is also 
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Fig. 3. Left Panel: Angular profiles of the expected DM induced synchrotron source 
(solid lines) and of the detected diffuse emissions (dotted lines) at 90 cm in the sur- 
veys of Refs. (23) (green) and (0) (black). The DM signal profile is shown also for 
a hypothetical EVLA observation with FWHM=200"(red). We consider as bench- 
mark case the A sp halo profile, Mdm = 3 TeV and a ann v = 3- 10~ 28 cm 3 s~ 1 . Dashed 
lines show the experimental sensitivities. Right Panel: For a few selected values of 
the DM mass, detectability of a monochromatic gamma-ray signature by the CTA 
project as a function of effective area x exposure time. The latter is expressed in 
terms of 1 km 2 x 50 hours, which can be considered as a conservative estimate for 
one year of observation by CTA. 

shown, together with the 3cr sensitivity of the detectors. We take as bench- 
mark case for the A_ candidate, a mass of 3 TeV and an annihilation rate 
of cr ann v = 3 • 10~ 28 cm 3 s _1 . The DM distribution considered is again the A sp 
profile. If the astrophysical radio diffuse emission is approximately isotropic 
at any scales, bounds on WIMP parameter space that could be extracted are 
not so stringent, as shown by the green curves, which is averaged over an 
angular resolution of 40 arcmin. On the other hand, if, on smaller scales, re- 
gions without contamination from astrophysical background are present, this 
type of surveys seems to be very promising, as shown in particular by the red 
curves, representing a hypothetical observation by EVLA with FWHM=200". 
However, this picture is probably based on a too optimistic assumption and 
it has to be considered limiting case. 

So far we have considered only continuum energy spectra of photons and elec- 
trons/positrons. The coupling between A_ and electrons is very tiny, since the 
latter are completely localized on the 4D brane at the boundary of the extra 
dimension. Thus, for our purposes, the prompt emission in monochromatic 
electrons and positrons can be neglected. Even if a gamma-ray continuum 
signal from DM annihilation exceeds the astrophysical background, the iden- 
tification of the DM component, which involves the exclusion of any other 
astrophysical explanation, could be a difficult task. The real "smoking gun" 
of a DM induced gamma-ray signal would be a monoenergetic spectral signa- 
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ture. By definition, the DM coupling with photons is highly constrained, but 
a direct WIMP annihilation into 77 at one-loop level is allowed, producing 
photons with energy ~ M DM , being WIMPs non-relativistic. Since A_ is 
an Abelian gauge boson, this process can occur through fermionic boxes. The 
main contribution is given by fermion triplets in the loop, for the same reason 
(i.e. the derealization) stated above referring to the tree level annihilation 
into fermions. The cross section computation is performed following Ref. |3ol ). 
and obtaining a 77 f ~ 2 • 10~ 4 <J ann v. The total number of events associated 
to DM annihilations into monochromatic 77 in a detector pointing to the GC 
direction with angular resolution Af2, is given by: 

-13 VlineV ( TeV \ 2 = A eff T 



N Hne = 1.9 10-" s m 3 -1 UfH ^(Afi)Afi^LL- . (2) 
10 61 cm a s 1 \M DM J m A s 

The quantity J(AQ), containing all the spatial information, is defined as: 

J(M]) = 8^c( o.3GeV/cm3 ) 2 M/^ eXp ( " 2^4) L/®* 

(3) 

where p is the DM halo profile, 9 is the angular off-set with respect to the 
GC and I is the coordinate along the line of sight. In Eq. [3] we consider a 
circular Gaussian angular resolution of width 9 a for the detector. The ratio 
between the gamma-ray signals originated in an A sp and an NFW profiles is 
given by the ratio: b = J^ sp (10~ 5 sr)/ J/vi?n/(10~ 5 sr), assuming AQ = 1CT 5 sr 
for modern ACTs. In the range of mass and cross section of the A_ model, it 
approximately follows the law: b ~ 10 4 [{a ann v /10 _28 cm 3 s _1 ) (TeV /Mdm)\ ° § - 
The dependence from the ratio a ann v / Mdm reflects the fact that the initial 
DM distribution, from which the A sp profile is derived, has a spike around 
the SMBH. In this case, self-annihilations frequently occur in the innermost 
region, triggering the final shape. 

The number of events associated to the 7-ray continuum background in a CTA 
bin can be obtained integrating the spectrum of the detected GC source and 
of the misidentified showers from hadrons and electrons (jsil) over an energy 
resolution of 10%. The probability of disentangling Ni events associated to the 
DM induced gamma-ray line from N^ g events of the continuum background is 
related to Odet = A"; / J N^ a + e 2 svs Nf , where e sys gives the level of systematic 
errors, taken to be 1% for CTA (I26l ). We estimate N[ to be a fraction €dm ~ 
2.7% of the total number of events. At fixed systematic error, the maximal 
significance which can be achieved increasing the effective area or the exposure 
time is a^ x = eoAi/tsys, i-e. the plateau in Fig. [3]b. A conservative guess for 
A e ffXT exp is 1 km 2 x 50 hours in one year of observation by CTA. As shown 
in Fig. [3]b, the prompt monochromatic emission of 77 originated from A_ 
annihilation in an A sp halo profile needs an extra factor of 100 in A e ffXT exp 
in order to be detected at ~ 3a; this could be reached only with a quite larger 
setup than the minimal designed and in several years of observation. 
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To conclude, in this Letter we have outlined the properties of a DM candidate 
recently proposed, sketching its multi-wavelength indirect signal from the GC. 
Cosmology and EW precision tests fix its mass and total annihilation cross 
section in a narrow window, which is compatible with the bounds associated 
to the detected emissions at the GC, but can be definitely tested by the forth- 
coming gamma-ray and wide-field radio surveys, if the Milky Way halo profile 
is spiky. We also discuss the possible detection of an induced gamma-ray line 
in the same framework. On the other hand, in case of NFW or more shal- 
low profiles, the model cannot be constrained through this multi-wavelength 
strategy. 

We would like to thank G. Bertone and D. Merritt for kindly providing some 
of DM halo profiles which have been used in this analysis. We also would like 
to thank M. Serone and P. Ullio for useful discussions. 
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